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METHODE  DE  CALCUL  DU  FLUTTER  EN  PRESENCE  DE  JEU  MECANIQUE 
ET  VERIFICATION  EXPERIMENT  ALE 

( FLUTTER  ANALYSIS  METHOD  IN  PRESENCE  OF  MECHANICAL  PLAY  AND  EXPERIMENTAL 

VERIFICATION) 

par 

C.  PETIAU,  B.  JOURNEE,  E.  GARRIGUES 
Aviation  Marcel  Dassault 
78,  Quai  Marcel  Dassault 
Cedex  300 

92552  St  Cloud  Cedex 
France 


RESUME 

Le  calcul  de  flutter  en  presence  de  jeu  mecanique 
intervient  pour  la  certification  des  timoneries  fail- 
safe des  empennages  des  Falcon. 

Une  methode  de  calcul  appropriee  a ete  developp^e 
dans  notre  logiciel  ELFINI ; elle  est  fondee  sur : 

- une  modelisation  structurale  dynamique 
Elements  Finis, 

- la  representation  des  forces  aerodynamiques 
instationnaires  dans  le  domaine  temps  par 
modele  type  Karpel, 

- une  integration  implicite  dans  le  temps, 

- la  resolution  des  jeux  et  contacts  par  un 
algorithme  derive  des  techniques  d'optimisation 
quadratique. 


Nous  exposons  la  methodologie  de  validation 
experimentale  sur  maquette  dynamique  en 
soufflerie,  en  3 etapes  : 

- definition,  calcul  preliminaire,  dimension- 
nement  de  l'experience  pour  verifier 
l'observabilite  des  phenomcnes  recherches, 

- verification,  calibration  du  modele  structural 
par  essais  statique  et  de  vibration,  du  module 
aerodynamique  k partir  de  mesures  de  pressions 
stationnaires  et  instationnaires  et  confrontation 
des  vitesses  critiques  de  flutter 
calculees/mesurces  contact  bloque, 

- confrontations  calculs/essais  de  flutter  en 
presence  de  jeux/contacts. 

En  procedant  a partir  des  modeles  structuraux  et 
aerodynamiques  recales  classiquement,  la 
simulation  restitue  les  resultats  d'essais  a la  fois 
qualitativement  (comportement  amorti,  cycles 
limites,  divergences  et  quantitativement  (niveaux 
d'acceleration). 


ABSTRACT 

Flutter  analysis  in  presence  of  mechanical  play 
occurs  in  the  certification  of  fail-safe  linkage  of 
Falcon  horizontal  stabilizer. 

The  corresponding  calculation  method  has  been 
developed  in  our  ELFINI  software,  it  is  based  on  : 

- a Finite  Element  structural  dynamics  model. 


a time  domain  model  of  unsteady  aerodynamics 
(Karpel  type). 


- an  implicit  time  integration, 

- the  resolution  of  play/contact  by  an  algorithm 
derived  from  quadratic  optimization  techniques 


We  expose  the  experimental  verification 
methodology  with  dynamic  model  in  wind  tunnel, 
in  3 steps : 

- definition,  preliminary  calculations,  sizing  of 
experience  for  checking  observability  of  studied 
phenomena, 

- verification,  calibration  of  structural  F.E.  model 
with  static  and  vibration  tests,  of  aerodynamic 
model  with  steady  and  unsteady  pressure 
measurements,  comparison  of  calculated  / 
measured  critical  flutter  speeds  with  fixed 
contact, 

- flutter  calculation/test  comparison  in  the 
presence  of  plays/contacts. 


From  classically  calibrated  structural  and 
aerodynamic  models,  the  simulation  reproduces 
test  results  both  qualitatively  (damped  behaviour, 
limit  cycles,  divergence)  and  quantitatively 
(acceleration  levels). 


Paper  presented  at  the  RTO  AVT  Specialists’  Meeting  on  "Structural  Aspects  of  Flexible  Aircraft  Control”, 
held  in  Ottawa,  Canada,  18-20  October  1999,  and  published  in  RTO  MP-36. 
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INTRODUCTION 

Le  probleme  de  l'analyse  du  flutter  en  presence  de 
jeu  mecanique  se  pose  regulierement  pour  la 
certification  des  avions  civils.  Dans  le  cas  des 
avions  DASSAULT  (Mercure,  Falcon),  il  resulte 
du  montage  "fail-safe"  de  la  timonerie  de  calage  de 
l'empennage  horizontal ; un  element  de  secours  est 
monte  avec  jeux  pour  ne  transmettre  les  charges 
qu'en  cas  de  defaillance  de  l’element  principal  (voir 
planche  1). 

Pour  traiter  ce  cas,  nous  avons  developpe  des  la  fin 
des  ann^es  70,  et  perfectionne  depuis,  un  module 
d'analyse  specifique  dans  le  cadre  de  notre  systeme 
d'analyse  structurale  et  aeroelastique  ELFINI.  Ce 
module  couple  l'analyse  dynamique  Elements  Finis 
avec  non  linearite  de  contact  avec  un  calcul 
d'aerodynamique  instationnaire  domaine  temps  ; 
nous  en  presentons  les  grandes  lignes  au  § 2. 


2.  PRINCIPES  DE  LA  METHODE  DE 
CALCUL  DE  FLUTTER  AVEC  JEU 

La  methode  a plusieurs  variantes  dont  les  points 
communs  sont : 

- une  modelisation  elastique  dynamique 
Element  Fini  classique  (voir  maillage  type 
planche  2) 

- une  reduction  de  base  (facultative,  voir  plus 
loin)  pour  limiter  le  nombre  de  degres  de 
liberte,  variantes  base  "modale"  ou  base  de 
"charge"  (principe  specifique  a ELFINI,  voir 
reference  1),  dans  tous  les  cas  la  base  doit  etre 
enrichie  des  deformees  statiques  pour  des 
charges  unitaires  aux  points  de  contact. 

Hors  modelisation  jeux/contact,  l'equation  de  la 
dynamique,  dans  la  base  de  calcul,  s'ecrit : 


La  certification  d'un  avion  fonde  sur  un  tel  outil  a 
pose  des  le  depart  le  probleme  de  sa  validation. 

On  a commence  par  des  verifications  partielles,  en 
s'assurant : 

- de  la  convergence  et  de  la  stabilite  des  solutions 
avec  les  hypotheses  de  modelisation 
(variantes/finesse  du  modele  E.F.,  nombre  et 
type  des  degres  de  liberte  dynamiques, 
raffinement  du  modele  aerodynamique), 


® [m]x’+[C]x*+  [k]x  = f.aero  + f.autres 

- Representation  rationnelle  des  forces 
a£rodynamiques  instationnaires  (approche 
type  Karpel,  ref.  2) 

Les  forces  aerodynamiques  instationnaires, 
calculees  normalement  par  ELFINI  dans  le 
domaine  frequence  sont  approximees  (lissage 
au  moindre  carre)  par  une  fonction  de  transfert 
de  la  forme  : 


- du  comportement  aux  limites  (contact 
permanent,  jamais  de  contact)  qui  doit  etre 
recoupe  par  les  calculs  classiques, 

- de  la  validite  des  composants  de  la  methode  de 
calcul,  qui,  etant  employes  dans  d'autres 
problemes,  ont  fait  par  ailleurs  l'objet  de 
verifications  experimentales  (exemple  contact 
statique,  impact  de  projectile,  interaction  avion 
atterrisseurs,  flutter  de  panneau  avec  non 
linearite  de  membrane  des  tuiles  de  la  navette 
Hermes, ...). 

Pour  une  verification  experimentale  "complete",  les 
essais  en  vol  etant  a l'evidence  trop  risques,  il  a ete 
decide  de  lancer  la  campagne  de  validation  par 
essais  sur  une  maquette  aeroelastique  dynamique 
en  soufflerie  dont  nous  presentons  la  demarche  et 
quelques  resultats  significatifs  (§  3). 

Du  fait  de  son  cout  important,  cette  operation, 
supportee  par  les  Services  Officiels  Fran^ais 
(SPAe),  a ete  menee  conjointement  entre 
DASSAULT,  AEROSPATIALE  et  l'ONERA. 


f.aero=[-i  pV‘A„+i  pVsA,+ipsJA2 


+|PVsD(^I-R)‘E]l 


correspondant  a une  representation  par 
equations  differentielles : 

® f.aero=I-pV2AoX+I-pVAix‘ 

+ fpA2  x“+fpVDy‘ 

— -Ry  = E x 

V J 

Les  matrices  Ai  representent  des  termes  de 
rigidite,  amortissement  et  masse  ajoutee,  les 
termes  en  y (eventuellement  negligeables) 
correspondent  a l'effet  des  degres  de  libertes 
"internes"  du  fluide,  pour  la  suite  nous 
convenons  de  rassembler  dans  x tous  les  degres 
de  liberty,  ce  qui  aboutit  en  reunissant  ® et  ©a 
la  modelisation  du  mouvement  par  un  systeme 
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d'equations  differentielles  ordinaires  du  2eme 
ordre  du  type  : 

® [M]x‘*  + [B]  x’  + [K]  x = f autres 

- Schema  d'integration  implicite  de  I'equation 
precedente  (Newmark,  Humbolt,  etc  ...). 

Les  termes  en  derives  sont  exprimes  comme 
combinaison  lineaire  de  l'etat  x(t)  inconnu  et 
des  etats  precedents  xt  - nAt,  I'equation 
differentielle  (3)  devient  un  systeme  d'equations 
lineaires  ordinaires. 

3ou4 

© [H]  X(t)  = f autres+  ]TBn  x(t  - nAt) 

i 

- La  resolution  des  jeux  et  contacts  est  traitee 
comme  solution  d'un  systeme  d'equations  et 
d'inequations  de  la  forme  : 

[H]  x = f + [N]  a 

Equations  d'equilibre  aeroelastique 

[N]t  X - r > 0 

Inequations  des  jeux 

a < 0 

Compression  aux  points  de  contacts 

Elle  est  mende  par  une  extension  de 
l'algorithme  dit  de  "gradient  spherique"  utilise 
dans  ELFINI  pour  la  rdsolution  des  contacts  en 
elasticity  lineaire. 

La  solution  est  recherchee  comme  combinaison 
lineaire  de  la  solution  sans  contact 

[H]'f 

et  des  deformees  pour  des  reactions  de  contact 
unitaire 

1H1'  INc] 

soit 

x = [H]1  [f  + [Nc]  a] 

Les  termes  a,  reactions  de  contact  (ou 
multiplicateurs  de  Lagrange),  correspondent 
aux  reactions  de  contact,  ils  doivent  etre 
negatifs  (compression)  ou  nuls. 

Ils  sont  obtenus  par  la  relation  de  jeu  nul  au 
contact 

[Nc]t  X-rc  = 0 soit 

a = [ [Nc]t  [H] 1 [Nc]]  1 [ [Nc]t  [H]  1 f-rc] 

II  est  a noter  que  la  non  symetrie  de  [H]  ne 
remet  pas  en  cause  la  validity  de  l'algorithme 


du  "gradient  spherique"  bien  que  ce  dernier  soit 
presente  habituellement  (d'oii  son  nom)  comme 
correspondant  a une  minimisation  sous 
contrainte  de  forme  quadratique  (rendue 
"spherique"  par  changement  de  repere). 

Autour  des  principes  de  base  qui  ont  ete 
presentes,  la  methode  a plusieurs  variantes  : 

- sur  la  modelisation  structurale  E.F.,  qui  peut 
etre  menee  sans  reduction  de  base,  en 
particulier  quand  on  est  en  presence  de  non 
linearites  structurales  type  grand  deplacement 
ou  plasticite  (voir  calcul  des  tuiles  de  la  Navette 
Hermes  ref.  3). 

- sur  la  modelisation  aerodynamique  et  ses 
possibility  de  recalage  experimental  (voir  ref. 
4),  par  le  couplage  direct  avec  calculs  C.F.D. 
(non  lineaires  ou  linearisees,  voir  ref.  5). 

3.  VALIDATION  DE  LA  METHODE, 
CONFRONTATION  EXPERIMENTALE 

L'organisation  de  cette  etude  a suivi  une  trame 
semblable  a celle  de  l'enchainement 
conception/calculs/essais  d'un  projet  d'avion,  voir 
planche  3. 

• lere  etane  : definition  de  1'experience  et 
de  la  maquette  aeroelastique,  calculs  preli- 
minaires  supportant  le  dimensionnement 
et  la  verification  de  1'observabilite  des 
phenomenes  rccherches. 

Apres  avoir  envisage  initialement  de  tester  un 
empennage  entier  (derive  + P.H.)  de  Falcon  10 
dans  la  soufflerie  SI  de  l'ONERA  de  Modane, 
on  s'est  limite  a ime  maquette  dynamique 
generique,  de  la  forme  d'un  empennage 
horizontal  d'avion  type  Airbus  dote  d'une 
gouveme  de  bord  de  fuite,  testee  dans  la 
soufflerie  transsonique  S2  de  l'ONERA  (voir 
planche  2). 

Le  caisson  du  plan  fixe  est  fixe  a la  paroi  par 
l'intermediaire  d'un  talon  situe  a l'emplanture 
du  caisson.  La  gouveme  est  guidee  en  rotation 
par  quatre  chamieres  fixees  sur  le  caisson  du 
plan  fixe.  Le  blocage  de  la  gouveme  en  rotation 
est  assure  par  l'intermediaire  d'une  lame  fixee 
en  extremite  de  l'axe  d'articulation.  Cette  lame 
est  bloquee  en  rotation  par  un  dispositif 
reglable  k butee  que  l'on  peut  deplacer  le  long 
de  la  lame.  Ce  dispositif  associe  a des  lames 
d'epaisseurs  difierentes  facilement  inter- 
changeables  permet  de  regler  la  frequence  de 
rotation  gouveme  sur  une  valeur  choisie.  On  a 
adjoint  a cette  maquette  une  palette  sur  le 
saumon  du  caisson  afin  de  pouvoir  introduire 
une  excitation  forcce.  Cette  palette  est 
commandee  par  un  verin  hydraulique  situe  dans 
le  caisson  du  plan  fixe. 
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L'instrumentation  d'essai  est  constitute  par  un 
ensemble  de  17  ponts  de  jauge,  de  20 
accelerometres,  de  5 fibres  optiques,  de  84 
capteurs  de  pression,  ainsi  que  de  deux 
potentiometres  pour  mesurer  la  rotation  de  la 
gouveme. 

Cette  maquette  possede  plusieurs  degres  de 
liberte  au  niveau  du  blocage  de  la  lame  en 
rotation.  Ce  blocage  est  realise  au  niveau  de  la 
lame  par  un  dispositif  mobile  le  long  de  la  lame 
qui  permet  de  modifier  la  longueur 
d'encastrement  de  la  lame.  II  existe  plusieurs 
lames  d'epaisseurs  differentes.  Ces  deux 
possibility  donnent  une  grande  latitude  pour  le 
reglage  de  la  frequence  de  rotation  gouveme 
qui  pilote  le  mtcanisme  de  flottement  ainsi  que 
la  Pi  critique. 

A ce  dispositif  de  reglage,  on  a ajoute  un 
dispositif  permettant  d'introduire  un  jeu  au 
niveau  de  l'encastrement  de  la  lame  pour 
etudier  les  comportements  non  lineaires  en 
presence  de  jeu  en  soufflerie.  Un  schema  de  ce 
dispositif  est  presente  planche  2. 

Le  predimensionnement  de  l'experience  a ete 
mene  avec  le  modele  Elements  Finis  presente 
planche  4. 

• 2eme  Etape  : calibration  des  modeles 

elastique  et  aErodynamique  classiques,  hors 
effet  de  jeux,  menee  en  sous-etapes. 

- Validation  du  modele  Elements  Finis 
statique  avec  les  correlations  calcul-essais 
des  reponses  de  jauges  de  contraintes  a 12 
chargements  statiques.  Cette  correlation  est 
relativement  satisfaisante  en  dehors  de  tout 
recalage  (voir  planche  5). 

- Validation  du  modele  dynamique  (modes 
propres).  Nous  n'avons  eu  a proceder  qu'a  la 
seule  calibration  du  coefficient 
d'encastrement  de  la  lame  sur  l'axe 
(independante  de  l'epaisseur  et  la  longueur 
de  la  lame)  pour  obtenir  une  confrontation 
calcul  essais  satisfaisante  (voir  planche  6). 

- Validation  du  modele  aerodynamique  (a 
Mach  0.6).  Elle  a ete  fondee  sur  la 
confrontation  du  calcul  aux  mesures  de 
pression  pour  des  excitations  par  braquage 
de  la  gouveme  en  stationnaire  et 
instationnaire  (planche  7).  Le  resultat  est 
conforme  a notre  experience  habituelle  de  la 
methode  des  doublets : 

. le  calcul  prevoit  correctement  les 
pressions  sur  le  caisson 
. les  pressions  de  bord  d'attaque  sont 
surestimees 


. les  pressions  sur  la  gouveme  et  le 
moment  de  chamiere  sont  surestimes 
d'un  facteur  2 

Cette  confrontation  a confirm^  le  facteur  de 
correction  de  0.5  des  pressions  gouvemes  qui 
avait  dejii  ete  introduit  "par  experience"  dans 
les  calculs  de  predimensionnement.  Cette 
correction  a ete  reconduite  pour  la  suite  des 
calculs. 

- Validation  du  modele  de  calcul  du  flutter 
hors  effets  de  jeux/contacts. 

Des  resultats  types  de  la  confrontation  calcul 
essais  sont  presentes  planche  8.  Les 
comparaisons  de  pressions  dynamiques  (Pi) 
critiques  et  des  Evolutions  des 
frequences/amortissements  sont  satisfaisantes. 

Nous  avons  aussi  mene  une  verification  par 
experience  numerique  du  modele  de 
rationalisation  des  forces  aerodynamiques 
instationnaires,  le  modele  semble  quasi  exact 
tant  pour  la  reconstitution  des  forces 
generalise  que  pour  la  vitesse  de  flutter  (voir 
planche  9). 

• 3eme  Etane  : confrontation  des  resultats  de 
flottement  non  linEaire. 

Nous  en  presentons  des  extraits,  planches  10  et 
11,  concernant  un  jeu  "centre"  de  0.1  mm 
(epaisseur  lame  4,5  et  1,5  mm  longueur 
130  mm). 

Le  calcul,  comme  les  essais,  permettent 
d'identifier  3 regimes  de  fonctionnement : 
stabilite,  cycles  limites,  divergences.  Dans  le 
cas  des  cycles  limites  les  niveaux  d'acceleration 
calcules/mesures  se  comparent  bien. 

A noter  que  le  domaine  de  pression  generatrice 
des  cycles  limites  ne  coincide  pas  avec 
l'intervalle  des  pressions  critiques  gouveme 
libre  - gouveme  bloquee  ; on  peut  rencontrer 
des  cycles  limites  en  dessous  de  la  Pi  critique 
gouveme  libre  et  des  cycles  limites 
"conditionnels"  au-dela  de  la  pression  critique 
gouveme  bloquee,  en  fonction  de  l'energie 
initiale  donnee  au  systeme. 

4.  CONCLUSIONS,  LECONS  A TIRER 

La  chaine  de  calculs  de  modelisation  du  flutter  en 
presence  de  jeu  mecanique  est  deja 
convenablement  validee  avec  cette  premiere 
confrontation  experimentale. 

Les  causes  d'imprecision  potentielle  des  calculs 
n'apparaissent  pas  devoir  venir  des  parties  les  plus 
originales  de  l'algorithme  comme  la  modelisation 
dynamique  des  contacts  ou  l'aerodynamique 
instationnaire  domaine  temps,  mais  plutot  des 
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parties  "classiques"  du  modele  intervenant  dejti 
dans  le  calcul  du  flutter  lineaire,  en  particulier : 

- le  modele  Elements  Finis  dynamique  qui  doit 
etre  verifie/recale  sur  les  essais  de  vibrations  au 
sol 

- le  modele  des  forces  adrodynamiques 
stationnaires  et  instationnaires  qui  doit  etre 
recale,  si  on  le  peut  sur  des  essais  en  soufflerie 
et  surtout  sur  les  essais  en  vol  (voir  ref.  4). 


C'est  d'ailleurs  sur  l'aerodynamique  que  nous 
attendons  les  principaux  progres  des  calculs  de 
flutter  lineaire  ou  non  lineaire  comme  ici  ; nous 
remplagons  progressivement  les  modeles  classiques 
de  potentiel  aerodynamique  traites  par  singularity 
par  des  methodes  Elements  Finis  "Euler" 
nonlineaires  ou  "linearisees"  (voir  ref.  5).  Nous 
envisageons  ainsi  une  prochaine  etape  de 
confrontation  calcul -cssai  dans  la  zone  du 
transsonique  avec  la  meme  maquette,  cela  pour  des 
configurations  sans  et  avec  jeu  mecanique. 
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PLANCHE 1 

ORIGINE  DE  L' ETUDE 


probleme  de  certification 
de  Fempennage 


Mercure 


- Falcon  900 


- Falcon  2000 


conception  fail-safe  de  la 
bielle  d’attaque  de 
I’empennage 

solution  technique 
comportant  un  jeu 

comportement  dynamique 
non-lineaire 


Dispositif 
de  commande 
non-lineaire 


axe  d’articulation 
principal 


Empennage 


PLANCHE  2 


MAQUETTE  AEROELASTIQUE  D’ EMPENNAGE  TYPE  ’’AIRBUS” 
ESSAYE  A LA  SOUFFLERIE  S2  DE  L’ONERA 


palette  d’excitation 


Structure  de  type  aeronautique 

dispositif  de  reglage  de  la 
frequence  de  rotation  gouverne 

palette  d’excitation 

encastrement  au  mur 

250  capteurs  ( pression,  jauge, 
accelerometre,  fibre  optique) 


caisson 


encastrement 
dispositif  de  securite 


gouverne 

axe 

gouverne 


BUTEE  A JEU 


ORGANISATION  DE  L’ETUDE 

SEMBLABLE  AU  PROCESSUS  CONCEPTION  / CALCUL  / ESSAIS  D'UN  AVION 
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VALIDATION  DU  MODELE  STATIQUE 
(MODELE  PRELIMINAIRE,  MODELE  FINAL) 


9-10 


VALIDATION  / RECALAGE  DU  MODELE  DYNAMIQUE 
(RECALAGE  DU  SEUL  COEFFICIENT  D'ENCASTREMENT  DE  LA  LAME) 
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MODELE  AERODYNAMIQUE  INSTATIONNAIRE 
COMPARAISON  ESSAIS  - CALCULS  DES  PRESSIONS 
(Calcul  integrant  les  effets  d'aeroelasticite  dynamique) 
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=>  CORRECTION  FACTEUR  0.5  SUR  LES  PRESSIONS  GOUVERNE 


LA  RATIONALISATION  DES  FORCES  AERODYNAMIQUES  INSTATIONNAIRES  EST  EXACTE 
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SYNTHESE  ETUDE  DE  FLUTTER  AVEC  JEUX 


Pi  critique  lineaire 


Amplitude  des  cycles  limites : 
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(jeu  0.1  mm  appui  a 130  mm) 


